The decay of homogeneous anisotropic turbulence 
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We present the results of a numerical investigation of three-dimensional decaying turbulence with 
statistically homogeneous and anisotropic initial conditions. We show that at large times, in the 
inertial range of scales: (i) isotropic velocity fluctuations decay self-similarly at an algebraic rate 
which can be obtained by dimensional arguments; (ii) the ratio of anisotropic to isotropic fluctuations 
of a given intensity falls off in time as a power law, with an exponent approximately independent 
of the strength of the fluctuation; (iii) the decay of anisotropic fluctuations is not self-similar, their 
statistics becoming more and more intermittent as time elapses. We also investigate the early stages 
of the decay. The different short-time behavior observed in two experiments differing by the phase 
organization of their initial conditions gives a new hunch on the degree of universality of small-scale 
turbulence statistics, i.e. its independence of the conditions at large scales. 



I. INTRODUCTION 



Decaying turbulence has attracted the attention of var- 
ious communities and is often considered in experimental, 
numerical and theoretical investigations [1-3] . It is in fact 
quite common that even experiments aimed at studying 
stationary properties of turbulence involve processes of 
decay. Important examples are provided by a turbulent 
flow behind a grid (see [4] and references therein) or the 
turbulent flow created at the sudden stop of a grid pe- 
riodically oscillating within a bounded box [5]. In the 
former case, turbulence is slowly decaying going farther 
and farther away from the grid and its characteristic scale 
becomes larger and larger (see [4] for a thorough exper- 
imental investigation). Whenever there is sufficient sep- 
aration between the grid-size Lj„ and the scale of the 
tunnel or the tank L$ 3> Li n , a series of interesting phe- 
nomenological predictions can be derived. For example, 
the decay of the two-point velocity correlation function, 
for both isotropic and anisotropic flows, can be obtained 
under the so-called self-preservation hypothesis (see [3] 
chapter XVI). That posits the existence of rescaling func- 
tions allowing to relate correlation functions at different 
spatial and temporal scales. By inserting this assump- 
tion into the equations of motion, asymptotic results can 
be obtained both for the final viscosity-dominated regime 
and for the intermediate asymptotics when nonlinear ef- 
fects still play an important role. 

The status of the self-preservation hypothesis and the 
properties of energy decay in unbounded flows are still 
controversial [2-4,6]. Systematic results on related prob- 
lems have been established recently, e.g. for non-linear 



models of Navier-Stokes equations as Burgers' equation, 
see e.g. [7], and for stochastic models of linear passive 
advection [8], both in unbounded [9-12] and bounded do- 
mains [13, 14]. 

Here, we investigate the decay of three-dimensional ho- 
mogeneous and anisotropic turbulence by direct numeri- 
cal simulations of the Navier-Stokes equations in a peri- 
odic box. Previous numerical studies have been limited 
to either homogeneous and isotropic turbulence [15, 16] 
or to shell models of the energy cascade [17]. 
The initial conditions are taken from the stationary en- 
semble of a turbulent flow forced by a strongly anisotropic 
input [18]. The correlation lengthscale of the initial veloc- 
ity field Li n is of the order of the size of the box Lq « Lj„. 
In the first part of this paper, we shall try to answer the 
following questions about the intermediate asymptotic 
regime of nonlinear decay: How do global quantities, 
such as single-point velocity and vorticity correlations, 
decay ? What is the effect of the outer boundary on the 
decay laws ? Do those quantities keep track of the initial 
anisotropy ? As for the statistics of velocity differences 
within the inertial range of scales, is there a recovery of 
isotropy at large times ? If so, do strong fluctuations get 
isotropic at a faster/slower rate with respect to those of 
average intensity ? Do isotropic and anisotropic fluctua- 
tions decay self-similarly ? If not, do strong fluctuations 
decay slower or faster than typical ones ? 
In the second part we study the early stages of the decay, 
with the aim of establishing a link between the small- 
scale velocity statistics in this phase and in the forced 
case. That will allow us to argue in favor of an "ex- 
ponents only" universality scenario, for forced hydrody- 
namic turbulence. 
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II. NUMERICAL SETUP 



B. Decaying runs 



A. The initial conditions 



The initial conditions are taken from the stationary 
ensemble of a forced random Kolmogorov flow [18]. For 
sake of completeness, we recall here some of the statisti- 
cal properties of this forced turbulent flow. We consider 
the solutions of the Navier-Stokes equations for an in- 
compressible velocity field v. 



dtv + (v ■ V)u = — Vp + vAv, 



(1) 



in a three-dimensional periodic domain. To maintain a 
statistically stationary state Eq. (1) had to be supple- 
mented by an input term / acting at large scales. This 
force was strongly anisotropic: / = (0,0, f z (x)) with 
fz(x) = Fx cos[2irx/L x + fafy)] + F 2 co$[4irx/L x + 2 (t)], 
constant amplitudes Fx, 2 and independent, uniformly 
distributed, <5-corrclated in time random phases (f>x,2(t). 
This choice ensured the statistical homogeneity of the 
forcing and thus of the velocity field. We simulated the 
forced random Kolmogorov flow at resolution 256 3 for 
time spans up to 70 eddy turnover times [18]. The vis- 
cous term was replaced by a second-order hyperviscous 
term — is A 2 v. We stored 40 statistically independent con- 
figurations that here serve as initial conditions for the 
decaying runs. 




(a) v z (tq) 



(b) ««(10td) 



(c) Uz(IOOto) 




(d) v x (tq) 



(e) v x (10t ) 



(f) v x (W0t ) 



FIG. 1. Two-dimensional sections of a typical velocity field 
at different times of decay. The three top images and the three 
bottoms are the grayscale plots of the velocity components 
parallel and transverse to the direction of the forcing, respec- 
tively. Note in the upper row, the presence of anisotropic 
structures which decay as time elapses. 



As turbulence decays, the effective Reynolds number 
Re = VrmsLxx/v decreases, while the viscous characteris- 
tic scale and time increase. To speed up the numerical 
time-marching, it is then convenient to use an adapta- 
tive scheme. We calculate periodically the smallest eddy- 
turnover time from the energy spectrum and set the time 
step as 1/100 thereof. The whole velocity-field configura- 
tion is then dumped for offline analysis at fixed multiples 
{0, 1, 10, 10 2 , 10 3 , 10 4 , 10 5 , 10 6 } r of the initial large-scale 
eddy turnover time tq = Lo/vj^. In Fig. 1 we show a 
two dimensional section in the plane x-z of the velocity 
components v z and v x - 



III. THE DECAY OF GLOBAL QUANTITIES 

A first hint on the restoration of isotropy at large times 
can be obtained by the two-dimensional snapshots in 
Fig. 1. After a few eddy turnover times, it is evident that 
large-scale fluctuations become more and more isotropic. 
To give a quantitative measure, we collect for each run 
the temporal behavior of the following one-point quanti- 
ties: 



E, 



v l (t)v j {t), 



Qij = Wi(t)wj(t). 



(2) 
(3) 



By 777 we denote the average over space coordinates only, 
whereas (■ • •) will indicate the average over both initial 
conditions and space. The symmetric matrices Eij(t) 
and Qij(t) are then diagonalized at each time-step and 
the eigenvalues Ex(t), E 2 (t), E 3 (t) and Qi(t), n 2 (t),n 3 (t) 
are extracted. Since the forcing points in a fixed direc- 
tion two eigenvalues are almost degenerate, say E 2 and 
E3, and strongly differ from the first one, E\. The typi- 
cal decay of Ei(t) and fij(t) for i = 1, ... ,3 is shown in 
Fig. 2. During the self-similar stage, t € [10, 10 6 ], the 
energy eigenvalues fall off as £{1,2.3} ~ t 2 , a s expected 
for the decay in a bounded domain [4, 15]. The enstrophy 
eigenvalues, ^{1.2,3} decay as i -12 / 5 . The dimensional ar- 
gument that captures these algebraic laws proceeds as fol- 
lows. The energy decay is obtained by the energy balance 
dE(t)/dt = —e(t), where we estimate E ~ v 2 ms (t) and 
e(t) - vf ms {t)/L - E 3 / 2 (t)/L , and obtain E(t) - i~ 2 
and e(t) ~ <~ 3 . As for the vorticity decay, we have 
f2(t) ~ (S v v) 2 /rj 2 , where 77 is the dissipative lengthscale 
and 8 v v is the typical velocity difference at separation 77. 
Assuming a Kolmogorov scaling S n v ~ e(rj) 1 / 3 r7 1 / 3 , and 
recalling that for a second-order hyperviscous dissipation 
e(t) - v(^v) 2 lff we obtain 77 - i 3 / 10 and S n v - i" 9 / 10 , 
whence ~ t" 12 / 5 . We now focus on the process of 
recovery of isotropy in terms of global quantities. Wc 
identify two set of observables 
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A tJ E(t) = 



Ei{t)-Ej(t) 
Ei(t)+Ei(t) 
fii(t) - %(t) 



(4) 



(5) 



IV. THE DECAY OF SMALL-SCALE 
FLUCTUATIONS 

A. The self-preservation hypothesis for anisotropic 
decay 



which vanish for isotropic statistics. Their rate of decay is 
therefore a direct measurement of the return to isotropy. 
The energy matrix E^ is particularly sensitive to the 
large scales while small-scale fluctuations are sampled by 
. .As it can be seen from Fig. 3, both large and small 
scales begin to isotropize after roughly one eddy turnover 
time and become fully isotropic (within statistical fluc- 
tuations) after 100 eddy turnover times. However, small 
scales show an overall degree of anisotropy much smaller 
than the large scales. 



^(1,2,3)0 










\ 

\ \ 








\ \ \ 

V \\ 






Ei(t) E (2i3} (t) 


\ -1 2/5 








\\\ 


















\ 


4 10" 2 10° 


10 2 


10 4 


10 6 




t 







FIG. 2. Log-log plot of the eigenvalues of energy and vor- 
ticity matrices vs. time, expressed in to unit. 
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FIG. 3. Log-log plot of the anisotropy content at the large 
scales (Ai2.E(i), Ai3_E(t), top curves) and the small scales 
(Ai2f2(t), Ai3fi(i), bottom curves) as a function of time, ex- 
pressed in ro unit. The large-scale (small-scale) anisotropy 
content is defined as the mismatch between the eigenvalues of 
the single-point velocity (vorticity) correlation. 



The observables which characterize the decay of small- 
scale velocity fluctuations are the longitudinal structure 
functions 

S (n) (r, t) = ([(« (x + r,t)-v (x, t)) • r]"} . (6) 

Those quantities depend not only on the modulus of the 
separation r but also on its orientation f , since the ve- 
locity statistics is anisotropic. A method to systemati- 
cally disentangle isotropic from anisotropic contributions 
in the structure functions is based on the irreducible rep- 
resentations of the SO (3) group [19]. In this approach, 
the observables (6) are expanded on the complete ba- 
sis of the eigenfunctions of the rotation operator. The 
SO(3) decomposition of scalar objects, such as structure 
functions, is obtained by projection on the spherical har- 
monics Yj m (f): 

oo j 

S(")(r,i)=£ E S$(r,t)Y jm (r). (7) 

j=0 m=-j 

Here, (r, t) denotes the projection of the n-th order 
structure function on the (j, m) SO (3) sector, with j and 
m labeling the total angular momentum and its projec- 
tion in the direction z, respectively. Another equivalent 
possibility is to look at the SO (3) decomposition of the 
PDF of the longitudinal velocity differences. In this case, 
denoting by V(A, r; t) the probability that the longitudi- 
nal incremement S r v = (v(r,t) — v(0,t)) ■ f be equal to 
A, we may project V(A, r; t) on the SO(3) basis: 

oo j 

T(A, r;t) = £ £ ?V A 5 *) W)- (8) 

3=0 m=-j 

The projection Vj m (r, A; t) plays the role of an "effective 
PDF" for each single SO (3) sector. (It should be how- 
ever remarked that only the isotropic probability den- 
sity Voo iXi A; t) has the property of being everywhere 
positive and normalized to unity with respect to the 
weight r 2 /{Ait).) Indeed, the projections of the longi- 
tudinal structure function on any sector (j, m) can be 
reconstructed from the corresponding V ]m i r , A; t) by av- 
eraging over all possible A's: 

S$(r,t) = J dAA n V jm (r,A;t) . (9) 

That establishes the equivalence between the decompo- 
sitions (7) and (8). 

The main points broached here are about the long-time 
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properties of the SO(3) projections, S ™^ (r,t). A simple 
isotropic generalization of the self-preservation hypothe- 
sis (see, e.g. Ref. [2]) amounts to writing: 

S$(r,t) = V<${t)f$(r/L jm (t)). (10) 

Note that vj^ (t) takes explicitly into account the fact 
that large-scale properties may depend in a nontrivial 
way on both (j, m) and the order n. Furthermore, Lj m (t) 
accounts for the possibility that the characteristic length 
scale depend on the SO(3) sector. 

In analogy with the observations made in the stationary 
case [18,20-26] we postulate a scaling behavior 

The time behavior is encoded in both the decay of the 
overall intensity, accounted by the prefactors a^™](i), and 
the variation of the integral scales Lj m (t). The repre- 
sentation (11) is the simplest one fitting the initial time 
statistics for t — and agreeing with the evolution given 
by the self preservation hypothesis in the isotropic case. 
The power law behavior for /-"^ (r/Lj m (t)) can be ex- 
pected only in a time-dependent inertial range of scales 
r](t) <C r <C L(t). As for the exponents appearing in 
(11), their values are expectedly the same as in the sta- 
tionary case. In the latter situation it has been shown 
that they are organized hierarchically according to their 
angular sector j [24] : 

c (n) < c (n) <... C (») <.... (12) 

Since the isotropic sector has the smallest exponent, at 
any given time and for a given intensity of the fluctuation 
(selected by the value of n) we have a recovery of isotropy 
going to smaller and smaller scales. Yet, deviations of the 
scaling exponents from their dimensional expectations 
make the recovery at small-scales much slower than what 
predicted by dimensional analysis [25,26]. Moreover, 
there are quantities which should vanish in an isotropic 
field and actually blow up as the scale decreases [25]. 
Concerning the time evolution, it seems difficult to dis- 
entangle the dependence due to the decay of a^™] (t) from 
the one due to the growth of the integral scale Lj m (t). 
Here, we note only that the existence of a running refer- 
ence scale, Lj m (t) introduces some non- trivial relations 
between the spatial anomalous scaling and the decaying 
time properties, and those relations might be subject to 
experimental verification. In our case, the fact that the 
initial condition has a characteristic lengthscale compa- 
rable with the box size, simplifies the matter. Indeed 
we expect that Lj m (t) w Lq, and the decay is due only 
to the fall off of the global intensity a^(t). Unfortu- 
nately, a shortcoming is that the width of the inertial 
range Lo/rjit) shrinks monotonically in time, thereby lim- 
iting the possibility of precise quantitative statements. 



B. Numerical results 



An overall view of the SO (3) projections at all resolved 
scales and for all measured decay times is presented in 
Fig. 4 for n = 2 and the isotropic, (j = 0, m = 0), sector. 
The same quantities are presented in Fig. 5 for the most 
intense anisotropic sector (j = 4, m = 0). 
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FIG. 4. Log-log plot of the isotropic component of the 
second order projection (r,t) vs. r for 7 decay times, 
t = 0, ro, 10r , 10 2 to, 10 3 to, 10 4 To, 10 s To (from top to bottom). 
The two straight lines correspond to the inertial range slope 
(r, i) ~ r 0,7 (top) and to the smooth differentiate slope 

Soo (r^) ~ r ' 2 (bottom). 
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FIG. 5. The same quantities as in the previous figure but 
for the anisotropic sector j — 4, m = 0, i.e. log-log plot of 
S$(r,t)v8.r. 

We notice that as time elapses the dissipative range 
erodes the inertial one, as a consequence of the growth 
of the Kolmogorov scale r){t). We notice in passing that 
the two curves corresponding to t = and t = t al- 
most coincide, i.e. even small scales are unchanged de- 
spite their typical eddy turnover times is much smaller 
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than To • This finding has some consequences that will be 
discussed at length in section V. A similar qualitative 
trend is displayed by the most intense anisotropic sector, 
(j = 4, m = 0) shown in Fig. 5, even though oscillations 
at small scales spoil significantly the scaling properties 
at small separations r. 
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FIG. 6. Log-log plot of S^l (r,t) vs. r, at t = To. Symbols 
refer to (j = 0, m = 0) (top) and to (j = 4, m = 0) (bottom). 



moments on the isotropic and an anisotropic sector at 
r = 80, within the inertial range. 
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FIG. 8. Log-log plot of the second-order moment 



S'2 ( r = 80, t) for the isotropic sector (j = 0, m = 0) and 



the anisotropic sector (j = 4,m = 0), vs. time, for a sepa- 
ration r in the inertial range. In the inset we plot the same 
curves for the fourth-order moment ( r = 80, t). 
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FIG. 7. The same quantities as in the previous figure but at 
the later time t = 10 2 to. The sectors shown are (j = 0, m = 0) 
and (j = 4, m = 0) (top to bottom), as in the previous figure. 

In order to assess the relative importance of isotropic 
and anisotropic contributions, we plot in Fig. 6 the SO(3) 
projections of some (j, m) sectors for the order n = 2 and 
a fixed time, t = tq. Fig. 7 shows the same quantities 
as in Fig. 6 but at a later time, t = 100r . Although 
the small-scale behavior of anisotropic sectors readily 
becomes rather noisy, the various contributions are or- 
ganized hierarchically and the isotropic contribution is 
dominant, as expected. 

Let us now analyze quantitatively the time-decay of the 
structure functions at a fixed separation. In Fig. 8 we 
show the long-time decay of the second and fourth-order 
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FIG. 9. Hierarchical organization of anisotropic fluctua- 
tions at long times. Log-log plot of the anisotropic projec- 
tions normalized by the corresponding isotropic projection 
(see text), at two fixed scales r = 80 and r = 40 (inset) 
for n — 2, 4, 6 in the anisotropic sector j = 4, m = 0. Symbols 

40 (star); II^ (empty box). 
~ 0.3. Same symbols in the 



read as follows : LT^ (full box); II 



The straight line is t 
inset. 



with x" 



We observe that the anisotropic sectors decay faster 
than the isotropic one, i.e. at a fixed scale there is a 
tendency towards the recovery of isotropy at large times. 
The relative rate of decay can be quantified by the ob- 
servable 



H (n) 



(r,t) 



5t } (r,t) 



t~ 



(13) 



In Fig. 9 we show IF-J (r, t) at r = 80 for structure 
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functions of order n — 2,4,6 and for (4,0), one of the 
most intense anisotropic sectors. In the inset we also 
plot the same quantities at the smaller scale r = 40. 
All anisotropic sectors, for all measured structure func- 
tions, decay faster than the isotropic one. The mea- 



sured slope in the decay is about x) 



(n) 



0.3 almost 



independent, within the statistical errors, of the order 
n. Note that these results agree with the simple pic- 
ture that the time-dependence in (11) is entirely carried 
by the prefactors a^(t) and the value of the integral 
scales Lj m (t) is saturated at the size of the box. In- 
deed, by assuming that large-scale fluctuations are al- 
most Gaussian we have that the leading time-dependency 



c (2n) • i (2) (2n-2) y-, . 

or CLj m is given by a^apo • * or the isotropic sec- 



tor, 4o n) 



(4 2 o h 



™, and plugging that in (13), we get: 

Ilg (r,t) ~ afl{t)/a { V{t) ~ i"*' with X \ ~ 0.3(±0.1) 
independent of n. The quality of our data is insufficient 
to detect possible residual effects due to Lj m (t), which 

would make Xj depend on n and j because of spatial 
intermittency. 

The interesting fact that we measure decay properties of 
the anisotropic sectors which are almost independent of 
the order of the structure functions indicate that we must 
expect some non-trivial time dependence in the shape of 
the PDF's Vj m (r, A; t) for j > 0. The most accurate 
way to probe the rescaling properties of Vj m (r, A; t) in 
time is to compute the generalized flatness: 



*fi(r,t) 



(Sfi{r,t)] 



(14) 



Were the PDF projection in the (j, m) sector self-similar 
for t To, then K^{r,t) would tend to costant val- 
ues. This is not the case for anisotropic fluctuations, as 
it is shown in Fig. (10). The curves K^(r,t) are col- 
lected for two fixed incrtial range separations, r = 80 
and r — 40 (inset), for two different orders, n = 4, 6 
and for both the isotropic and one of the most intense 
anisotropic sectors (j = 4, m = 0) . The isotropic flatness 
tends toward a constant value for large t. Conversely, 
its anisotropic counterparts are monotonically increasing 
with t, indicating a tendency for the anisotropic fluctu- 
ations to become more and more intermittent as time 
elapses. A consequence of the monotonic increase of in- 
termittency for large times is the impossibility to find a 
rescaling function, g(t, r), which makes the rescaled PDF 

g(t,r)Vj m (r, g £ r ) i t) time-independent at large times. 
Let us notice that also the behavior in Fig. 10 is in qual- 
itative agreement with the observation previously made 
that all time dependencies can be accounted by the pref- 
actors ag(i). Indeed, assuming that the length scales 
Lj m (t) have saturated and that the large scale PDF is 
close to Gaussian, it is easy to work out the prediction 
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FIG. 10. Log-log plot of the generalized flatness, (r, t) 
of order n = 4, 6 for both the isotropic (two bottom curves), 
and the anisotropic sector (j = 4, m = 0) (two top curves) 
at r = 80, and as a function of time. In the inset we plot 
the same quantities, in the same order, at a different inertial 
range scale, r = 40. 



= X*(n/2-l). 



We conclude this section by a brief summary of the re- 
sults. We have found that isotropic fluctuations persist 
longer than anisotropic ones, i.e. there is a time-recovery, 
albeit slower than predicted by dimensional arguments, 
of isotropy during the decay process. We have also found 
that isotropic fluctuations decay in an almost self-similar 
way whilst the anisotropic ones become more and more 
intermittent. Qualitatively, velocity configurations get 
more isotropic but anisotropic fluctuations become, in 
relative terms, more "spiky" than the isotropic ones as 
time elapses. 



V. SHORT-TIME DECAY 



Let us now move to the properties of decay at short 
times (t <C t ). Universality of small-scale forced tur- 
bulence is at the forefront of both theoretical and ex- 
perimental investigation of real turbulent flows [2]. The 
problem is to identify those statistical properties which 
are robust against changes of the large-scale physics, that 
is against changes in the boundary conditions and the 
forcing mechanisms. Our goal here is to relate the small- 
scale universal properties of forced turbulent statistics 
to those of short-time decay for an ensemble of initial 
configurations. An immediate remark is that one cannot 
expect an universal behaviour for all statistical observ- 
ables as the very existence of anomalous scaling is the 
signature of the memory of the boundaries and/or the 
external forcing throughout all the scales. Indeed, the 
main message we want to convey here is that only the 
scaling of both isotropic and anisotropic small-scale fluc- 
tuations is universal, at least for forcings concentrated 
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at large scales. The prefactors are not expected to be 
so. There is therefore no reason to expect that quantities 
such as the skewness, the kurtosis and in fact the whole 
PDF of velocity increments or gradients be universal. 
This is the same behavior as for the passive transport of 
scalar and vector fields (see [8] and references therein). 
For those systems both the existence and the origin of the 
observed anomalous scaling laws have been understood 
and even calculated analytically for some instances in 
the special class of Kraichnan flows [27] . Here, it is worth 
stressing that the universal character of scaling exponents 
is shared by both isotropic and anisotropic fluctuations 
[28]. 

For the Navier-Stokes case we have a huge amount of 
experimental and numerical indications that the velocity 
field shows anomalous scaling [2], suggesting the exis- 
tence of phenomena "similar" to those of the linear case. 
However, carrying over the analytical knowledge devel- 
oped for linear hydrodynamical problems involve some 
nontrivial, yet missing, steps. For the Navier-Stokes dy- 
namics, linear equations of motion surface again but at 
the functional level of the whole set of correlation func- 
tions. In a schematic form: 

d t C {n) = r(" +1 )c(" +1 ) + vD {n) C (n) + F (n \ (15) 

where r(" +1 ) is the intcgro-differential linear operator 
coming from the inertial and pressure terms and C^ n+1 ^ 
is a shorthand notation for a generic (n+l)-point correla- 
tor. The molecular viscosity is denoted by v and is 
the linear operator describing dissipative effects. Finally, 
is the correlator involving increments of the large- 
scale forcing / and of the velocity field. Balancing in- 
ertial and injection terms gives dimensional scaling, and 
anomalously scaling terms must therefore have a different 
source. A natural possibility is that a mechanism similar 
to the one identified in linear transport problems be at 
work in the Navier-Stokes case as well. The anomalous 
contributions to the correlators would then be associated 
to statistically stationary solutions of the unforced equa- 
tions (15). The scaling exponents would a fortiori be 
independent of the forcing and thus universal. As for the 
prefactors, the anomalous scaling exponents are positive 
and thus the anomalous contributions grow at infinity. 
They should then be matched at the large scales with 
the contributions coming from the forcing to ensure that 
the resulting combination vanish at infinity, as required 
for correlation functions. Our aim here is not to prove 
the previous points but rather to check over the most ob- 
vious catch: the Navier-Stokes equations being integro- 
differential, non-local contributions might directly couple 
inertial and injection scales and spoil the argument. This 
effect might be particularly relevant for anisotropic fluc- 
tuations where infrared divergences may appear in the 
pressure integrals [29]. 



1.5 




t 

FIG. 11. Top: Temporal decay of the second-order 
isotropic structure function S^(r, t), rescaled by its value at 
t = 0. Here r = 30, inside the inertial range. The two curves 
refer to the time evolution of the structure function starting 
from the forced-stationary velocity fields (squares, set A) and 
from the randomly dephased velocity fields (circles, set B). 
Time is normalized by the integral eddy turnover time. No- 
tice that for set B we observe changes on a time scale faster 
than the integral eddy turnover time. That is to be contrasted 
with the case A, where structure functions are strictly con- 
stant in time up to an integral eddy turnover time. Bottom: 
The same curves but for the fourth-order structure function. 

In order to investigate the previous point, we per- 
formed two sets of numerical experiments in decay. The 
first set, A, is of the same kind as in the previous section, 
i.e. we integrated the unforced Navier-Stokes equations 
(1) with initial conditions picked from an ensemble ob- 
tained from a forced anisotropic stationary run. Statisti- 
cal observables are measured as an ensemble average over 
the different initial conditions, (•) ens . The ensemble at 
the initial time of the decay process is therefore coincid- 
ing with that at the stationary state in forced runs. If cor- 
relation functions arc indeed dominated at small scales 
by statistically stationary solutions of the unforced equa- 
tions then the field should not decay. Specifically, the 
field should not vary for times smaller than the large-scale 
eddy turnover time r ~ L / '{v 2 ) 1 / 2 , with L denoting 
the integral scale of the flow. Those are the times when 
the effects of the forcing terms start to be felt. Note that 
this should hold at all scales, including the small ones 
whose turnover times are much faster than t . The sec- 
ond set of numerical simulations (set B) is meant to pro- 
vide for a stringent test of comparison. The initial con- 
ditions are the same as before but for the random scram- 
bling of the phases: Vi(k) — ► Pu(k)vi(k) exp(i#;(fc)). 
Here, v denotes the Fourier transform and Pu(k) is the 
incompressibility projector. In this way, the spectrum 
and its scaling are preserved but the wrong organization 
of the phases is expected to spoil the statistical stationar- 
ity of the initial ensemble. As a consequence, two differ- 
ent decays are expected for the two sets of experiments. 
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In particular, contrary to set A, set B should vary at 
small scales on times of the order of the eddy turnover 
times T r ~ r 2 / 3 . This is exactly what we found in the nu- 
merical simulations, as can be seen in Fig. 11, where the 
temporal behavior of longitudinal structure functions of 
order 2 and 4 is shown. The scaling of the contributions 
responsible for the observed behaviour at small scales are 
thus forcing independent. 

As for anisotropic fluctuations, we also found two very 
different behaviors depending on the set of initial con- 
ditions. In Fig. 12 we show the case of the projection 
Sjm( r — 60, i) for the anisotropic sector j — 4,m = 0. 
As it can be seen, for set A of initial conditions the func- 
tion is indeed not decaying up to a time of the order of 

TO- 

To conclude, the data presented here support the conclu- 
sion that nonlocal effects peculiar to the Navier-Stokcs 
dynamics do not spoil arguments on universality based 
on analogies with passive turbulent transport. The pic- 
ture of the anomalous contributions to the correlation 
functions having universal scaling exponents and non- 
universal prefactors follows. 



situations where the integral scale cannot grow indefi- 
nitely. The decay process at long times is then governed 
by the set-up at large scales. Anisotropies decay in time 
at a rate almost independent of the order of the moments 
and of the kind of anisotropic fluctuations. Projections 
of the PDF's on different SO (3) sectors show different 
intermittent properties (with the anisotropic sectors be- 
ing more intermittent). An obvious further development 
of this study would be to investigate the case where the 
integral scales Lj m (t) vary in time. Additional intcrmit- 
tency in time might then be brought in by the anomalous 
scaling in the space variables of the correlation functions. 
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FIG. 12. The same curves as in the previous figure but for 
the anisotropic sector (j =4, m = 0). 



VI. CONCLUSIONS 

We have presented a numerical investigation of the de- 
cay of three-dimensional turbulence, both isotropic and 
anisotropic. Concerning short-time decay, we have com- 
pared the decay for two different sets of initial conditions, 
with and without phase correlations. That gave some 
new hints on the properties of universality of isotropic 
and anisotropic forced turbulence. As for long times, 
we have found that fluctuations in the inertial range 
become more and more isotropic. On the other hand, 
the anisotropic components become more and more in- 
termittent, i.e. relatively intense anisotropic fluctuations 
become more and more probable. The main issue here 
was to investigate the effects of a bounded domain, i.e. 
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